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QCD strings originate from high-energy scattering in the form of Reggeons and Pomerons, and 
have been studied in some detail in lattice numerical simulations. Production of multiple strings, 
with their subsequent breaking, is now a mainstream model of high energy pp and pA collisions. 
Recent LHC experiments revealed that high multiplicity end of such collisions show interesting col¬ 
lective effects. This ignited an interest in the interaction of QCD strings and multi-string dynamics. 
Holographic models, collectively known as AdS/QCD, developed in the last decade, describe both 
hadronic spectroscopy and basic thermodynamics, but so far no studies of the QCD strings have 
been done in this context. The subject of this paper is to do this. First, we study in more detail the 
scalar sector of hadronic spectroscopy, identifying “glueballs” and “scalar mesons,” and calculate 
the degree of their mixing. The QCD strings, holographic images of the fundamental strings, thus 
have a “gluonic core” and a “sigma cloud.” The latter generates a exchanges and collectivization 
of the strings, affecting, at a certain density, the chiral condensate and even the minimum of the 
effective string potential, responsible for the very existence of the QCD strings. Finally, we run 
dynamical simulations of the multi-string systems, in the “spaghetti” setting approximating central 
pA collisions, and specify conditions for their collectivization into a black hole, or the dual QGP 
hreball. 

PACS numbers: 11.25.Tq, 25.75.-q, 12.38.Mh 


I. INTRODUCTION 
A. QCD Strings and Holographic Models 

High energy hadronic collisions in 1960’s led to discov¬ 
eries of Regge phenomenology of scattering amplitudes 
and striking Regge trajectories connecting mesonic and 
baryonic states. The leading singularity, producing a 
term in the amplitude 

A(s,t) ~ (1) 

is known as the Pomeron, named after I. Pomeranchuk. 
The Pomeron dominates the high-energy cross sections 
and elastic amplitudes. 

These observations were then explained in terms of 
the QCD strings. In particular, the slopes of the Regge 
trajectories, including that of the Pomeron, a', were re¬ 
lated to the string tension. (Creators of the fundamen¬ 
tal string theory have kept this notation for the funda¬ 
mental string scale). String tension is also the slope of 
the linear potential of confinement, which has been stud¬ 
ied in detail in lattice numerical simulations. Its value, 
Ts = (420 MeV)^, is in fact used as a standard definition 
of “physical units” in various confining theories. 

QCD strings have been with us for about half century, 
and yet, the interest to them was rather unsteady. From 
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the lattice we learned that the QCD strings are surpris¬ 
ingly thin, with an r.m.s. width of only about 0.17 fm. 
The dynamics of these strings are well described by the 
simplest Nambu-Goto action. The exponential growth of 
density of states of strings produces the Hagedorn phe¬ 
nomenon, a rapid excitation of the string, responsible for 
the deconfinement transition. 

The discovery of QCD gave rise in 1970’s to weak cou¬ 
pling, or pQCD, methods. In the field of hadronic colli¬ 
sions, the “hard,” or BFKL, Pomeron has been derived 
through the re-summation of gluonic ladders. This ap¬ 
proach focused on the Pomeron intercept a(0), while ig¬ 
noring a “stringy” Q!'(0). 

We return to the recent derivation of the Pomeron 
below, and now switch to strong coupling models, 
which came into existence after seminal discovery of the 
AdS/CFT duality. This relates the 4-dimensional, A/” = 4 
super Yang-Mills theory at strong coupling and large 
number of colors with the weakly coupled Type IIB su¬ 
pergravity in an AdS^ x S'® background created by a set 
of Ds branes, [1]. 

Maldacena in [1] calculated the potential of two static 
charges connected by a string by deriving the shape of 
a bending (geodesic) string in the bulk spacetime. As 
the setting is scale invariant, the modified strong cou¬ 
pling Coulomb law remains V ~ l/r, with only the coef¬ 
ficient modified. The calculation for the charges mov¬ 
ing with fixed velocity ±u away from each other has 
been performed by Lin and Shuryak [2], who found that 
for small enough v there is a scaling solution for the 
falling string, generalizing Maldacena case into a “gen¬ 
eralized Ampere law.” Yet, above certain critical veloc¬ 
ity this becomes unstable, and the stable solution in 
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that range has been found numerically. The hologram 
of that string calculated in [3] showed a near-spherical, 
non-hydrodynamical explosion, demonstrating quite di¬ 
rectly that in A/'=4 at strong coupling there should be no 
jets, even in e“'"e“annihilation into quarks. Consequently, 
suppose one can describe a set of falling strings in AdS^ 
(a “strongly coupled glasma”). If the endpoints are sep¬ 
arated by a very large rapidity interval (say ~ A?/=15 at 
the LHC), one would expect for the falling strings simi¬ 
larly large rapidities of the transverse flow. 

Both the gauge theory in question and the bulk AdS^ 
spacetime possess conformal symmetry. While vital for 
the discovery of the duality, conformality is certainly 
not welcomed when a description of the real world is 
intended. Over the years, the AdS/QCD approach has 
been developed, which includes violation of conformality. 

Unlike lattice QCD, AdS/QCD is not restricted to the 
Euclidean domain, and thus various time-dependent pro¬ 
cesses can also be studied. This opens a door to studies 
of various out-of-equilibrium settings devoted to under¬ 
standing of the matter equilibration. For example, mat¬ 
ter equilibration can be described via a “shock wave in 
the 5**' dimension”, which is equivalent to a membrane 
falling under its own weight [4] . 

So far, the AdS/QCD approach has been used to model 
the hadronic spectrum, as well as thermodynamics of 
QCD at finite temperatures. A particularly successful 
model of this class is the Improved Holographic QCD 
model, which incorporates several low energy QCD fea¬ 
tures like confinement and linearity of glueball trajec¬ 
tories, and reproduces the thermodynamic functions at 
finite temperature [5-8]. The flavor degrees of freedom 
are usually studied in the probe approximation, where 
they do not backreact to the glue [9-12]. Important re¬ 
cent advances [13], which are used below, include the 
back-reaction of the quarks to the glue in the Veneziano 
limit of QCD (also called V-QCD), in which the number 
of colors and flavors are comparable: 

Nc, Nf ^ oo Nf/Nc = X = const. (2) 

Here, we are interested in the QCD-string interaction 
and the effect of the QCD strings on the chiral symme¬ 
try where cr-meson plays a crucial role since it mainly 
mediates the string interaction and it is responsible for 
chiral symmetry restoration, as explained in (HI A, IV). 
Hence, it is necessary to work in the Veneziano limit 
where the flavor degrees of freedom are not suppressed. 
In the holographic context, the string directly sources 
the scalar glueballs, which are coupled to the cr-meson to 
leading order in the Veneziano limit [14, 15]. 

B. Strings in High Energy Collisions 

As has been already mentioned, high-energy collisions 
of hadrons are dominated by the Pomeron physics. For 
recent derivation of the Pomeron amplitude in string- 
based holographic settings, see [16, 17]. Without going 


into details, we remind the reader that the Pomeron pro¬ 
cess is basically a spontaneous reconnection of the strings; 
the hadrons, instead of being color neutral, are connected 
by (at least) two QCD strings after the collision. At the 
moment of the collisions, the strings are located in the 
transverse plane, appearing from a tunneling - under the 
barrier - process [16]. 

High energy colliders - RHIC at BNL and LHC at 
CERN - continue to provide new data on a wide range 
of hadronic collisions, from the basic proton-proton (pp), 
to proton-nuclei {dAu and pAu, to be called pA below), 
and heavy-ion AA collisions. As the multiplicity grows 
in each of the systems, a transition is observed. In the 
AA case, peripheral collisions are superpositions of NN 
collisions, while more central collisions lead to production 
of the QGP fireball, which subsequently explodes, leading 
to a set of observed collective phenomena, such as radial, 
elliptic, etc. flows. The typical pp, pA events are, on 
the contrary, well explained by pQCD and QCD strings, 
decaying independently, according to Lund-type models, 
such as Pythia and its descendants. 

The question where the transition between those two 
regimes happens is currently under intensive study, and 
collective phenomena such as radial, elliptic, and trian¬ 
gular flows were indeed observed in high multiplicity pA, 
and perhaps even pp collisions. 

For recent review of this phenomenology, one can see 
e.g. [18]. Let us mention the main details of pA collisions. 
Note that the typical impact parameter, 

6~y^^«1.6fm, (3) 

(in the LHC energy range) is about one order of mag¬ 
nitude larger than the string width of .17 fm extracted 
from the lattice data. The diluteness of the QCD strings 
in such a state is given by the ratio 

d = Ns(^'^^y ^10-^N,. (4) 

So, we see that even for Ng ~ 40 in the central pA colli¬ 
sion at the LHC, d < 1 and the strings are not yet fully 
overlapping. 

The central issue - referred to as the “radial flow 
puzzle” - explains why we are so interested in pA col¬ 
lisions. The initial density (multiplicity per transverse 
area dN/dAj_) suggests the following order of magnitude 
of various processes, 

dN^Xheral dN^^tral . dK^hest 
dAj_ (iAj_ dAi_ dAj_ ’ 

( 5 ) 

and, naively, one expects that the magnitude of the ex¬ 
plosion (the velocity or rapidity v± = tanh(i/_L) of the 
radial flow) would also follow this pattern. And yet, ex¬ 
periment gives us a different trend, 

AA,peripheral ^ AA,central ^ pA,central pp,hiqhest 

2/_l < 2/x < 2/x ~ V ± 

( 6 ) 
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Central pA collisions are therefore more “explosive” than 
expected. 

One scenario proposed to explain this puzzle, proposed 
by Kalaydzhyan and Shuryak [19], is a collective collapse 
of the string system created in maximal multiplicity pA 
case. If so, the size of the system is reduced and density 
is increased by a significant factor, leading to a different 
order of the densities, corresponding to the strength of 
the flow. Another notable consequence of this collapse 
is that the combined field of all of the strings becomes 
strong enough to restore chiral symmetry and thus create 
a QGP fireball, needed for an explosion. 

As time goes on and the remnants of the nucleons 
move forward and backward with large rapidities, the 
ends of the strings remain attached to them. As a result, 
the strings get extended longitudinally, and eventually 
break via quark pair production. The details of this pro¬ 
cess are modeled by descendants of the Lund model of 
particle production, such as the Pythia event generator. 
Since quark pair production proceeds via tunneling in 
the Schwinger mechanism, it is numerically suppressed. 
This allows a certain proper time Threaking 1-5 fm for 
strings to be stretched along the beam axis, forming a 
multi-string state. 

Since the creation of QGP fireball is dual to formation 
of a black hole in a bulk, this topic has been addressed 
in some of these works. We, however, focus instead on 
smaller systems such as typical pp or pA collisions, and 
thus instead of colliding walls, we discuss a string setting. 
Looking for a more amenable geometry, we look at the 
“spaghetti” phase, in which certain number of parallel 
QGD strings are produced and subsequently decay. The 
smallest number, originating from a single Pomeron or 
color exchange, is 2 strings, connected to leading quarks 
(diquarks). 

It has been argued in [19] that when the number of 
strings in the “spaghetti” gets large enough, mutual at¬ 
traction between them grows enough to induce a collapse 
and the collectivization of their fields. It has been argued 
that it should happen for the largest number in observed 
high multiplicity pA events, in which perhaps 40 or so 
strings are produced. In this paper, we extend this sce¬ 
nario into a holographic setting, where we calculate the 
string potential from the underlying holographic model. 


II. THE SETTING 
A. Background Solutions 

The IHQGD model was developed by Kiritsis et al. 
in [5-8]. This is a phenomenological holographic model 
of 4-dimensional Yang Mills at large N^- It includes an 
asymptotically AdS^ metric with logarithmic corrections 
corresponding to the perturbative running of the’t Hooft 
coupling at high energies. The model also contains a 
scalar dilaton field, which is the holographic analogue of 
the’t Hooft coupling. A sophisticated choice of the dila¬ 


ton potential leads to the correct description of certain 
low-energy features of the theory. The mesonic sector of 
QGD is modeled by placing Nf coincident pairs of II4-l)4 
branes in the background metric. As was shown in [ 11 ] 
and further analyzed in [ 12 ], chiral symmetry is broken 
as soon as the background metric is confining. 

The IHQGD model was generalized in [13] by includ¬ 
ing the back-reaction of the quarks, in a model named 
V-QGD. Its thermodynamics were studied in [20] and 
different aspects of its spectra in [14, 15, 21, 22]. For 
convenience of the reader, the definitions and terminol¬ 
ogy of the fields involved are explained in Table 1. 

The action for the metric and the dilaton $ is 

Sg = M^Nl j + (7) 

where M is the 5-dimensional Planck mass. The overall 
setting includes a background with a metric of the form 

= eyip{2A{z))[dz^ + gijdx^'dx^] , (8) 

where pij = diag(—,-I-,-I-,-f) is the Minkowski metric. 
The ’t Hooft A coupling is directly related to the dila¬ 
ton: A = exp($). The main field for the description of 
the flavor part of the theory is the tachyon (t), which is 
a bifundamental scalar under the SU{Nf)L x SU{Nf)]i 
flavor symmetry and is dual to qq operator. The back¬ 
ground fields are shown in Fig. 1. 

The action describing the vacuum of the flavor sector 
has been analyzed in [13] and is of the form 

Sf = -M^NcNf j d^xVf{X, t)^J det{g^i, + K,{\)d^Td^T ), 

(9) 

where and k,{X) are potentials which were con¬ 

strained by matching to QGD features in [15]. Here we 
use the class of potentials I, which are defined in Eqs. 
(4.19) and (4.20) of [15], with 17(0,0) = 3/11. 

The UV expansions of the background fields are such 
that they reproduce the perturbative running of the ’t 
Hooft coupling. 

log j + log(zA[/y) ’ 

t{z) ^ imqz{- log(A; 7 y z)y +£<7q z^{- log(Ac/y z))~'^ , 

( 10 ) 

where i is the AdS radius and Aj/y is the characteristic 
scale of the theory that corresponds to the QGD strong 
coupling scale. The parameter 7 is determined by the 
anomalous dimension of qq, [13], vriq is the bare quark 
mass, and Uq is the condensate. In the present work, we 
consider niq = 0. The determination of Uq is numerically 
very challenging; see [22]. The IR asymptotics of the 
fields are 

A{z) ~ —Aj^z^-I--log(A/flz) , ^{z) ~ , (11) 

T(z)^Toe^^^«^ ( 12 ) 
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FIG. 1: The background solutions for the metric scale factor 
(solid), the dilaton (dashed) and the tachyon (dotted) as a 
function of the holographic coordinate 2 . 


where the IR scale A/^ is an integration constant of the 
equations of motion. 

The scalar fluctuations of the model describe the mixed 
scalar glueballs and mesons and were analyzed in [15]. 
The relevant excitations of the tachyon, dilaton, and the 
scalar part of the the gij component of the metric are 

= $ + T = t + s, Qij = + . (13) 


B. Strings 

The bulk theory is argued to be a low energy effec¬ 
tive model of string theory, even if we are not in a con¬ 
trol string theory limit, so there are fundamental strings 
in the bulk. The electric fluxes running through these 
strings are sourced by charges at the boundary - the 
quarks of the gauge theory. In holographic models, stud¬ 
ies of stretching strings were first done in the original 
AdS/CFT correspondence with a conformal AdS^ space 
in refs. [2, 23]. In this setting, the middle section of 
the string is falling under gravity away from the bound¬ 
ary. The calculated holographic image of such stretching 
strings reveals a non-hydrodynamical explosion, which 
does not have any “jet-like” features (forward-backward 
peaks). 

The bulk strings in AdS/QCD background are not in¬ 
definitely falling, but can reach certain equilibrium posi¬ 
tions and “levitate”. Their hologram at the space bound¬ 
ary is what one would call the QCD strings. Thus in such 
models the forces between quarks gradually change from 
Coulomb 1/r at small r to linear R ~ r at large r, cor¬ 
responding to confinement, [6]. 

String dynamics are governed by the Nambu-Goto ac¬ 
tion 


Their invariant combinations under the linearized 5- 
dimensional diffeomorphism symmetry of the bulk space- 
time which correspond to the physical scalars are 

A' A' 

c = -^x , ^ = '^-^S. (14) 

C and ^ are dual to the RG operators of the boundary held 
theory which generate the glueball and scalar mesons, 
respectively, as a: —)■ 0. In the Veneziano limit (x = 
hnite), the two excitations mix to leading order, hence 
the distinction to glueballs and mesons requires a detailed 
study of the mixed fluctuation equations, (26, 27). We 
perform this analysis in section (IIIB) for x = 1 and 
make contact with phenomenology. We show that the 
mixing of the different huctuations can be weak even for 
hnite value of x. 


TABLE I: The fields and fluctuations of our model. 


4> = log A 

Dilaton 

A 

Metric conformal factor 

0 

+ 

II 

String frame metric conformal factor 

T 

Tachyon {qq scalar) 


Dilaton fluctuation 

s = St 

Tachyon fluctuation 

^ 

Scalar part of metric fluctuation 

C = i’- wx 

Scalar glueballs as a: —> 0 


Scalar mesons as a; —>■ 0 


Sng =-Tf j dTd<7^/-Aeigs , (15) 

igs)a/3 = {gs)iivdaX^dpX'' , (16) 

where we use the string frame metric 

( 5 s)/n/ = , As{z) = A{z) . (17) 

In top-down holographic models, the fundamental 
string tension Tf provides the input scale for the whole 
construction and is determined by the original string the¬ 
ory. Unfortunately, in the V-QGD model, such a connec¬ 
tion of Ty to other parameters is missing, and thus should 
be htted to phenomenology. 

The force of gravity - gradients of the metric - causes 
all objects, including strings, to fall to the IR (large z). 
The dilaton gradient, however, produces the opposite ef¬ 
fect. Specihcally, in the background used, the metric at 
large 2 is decreasing as A ~ —Airz^, but the $ contri¬ 
bution cancels this term; see Eq. (II). The conformal 
factor of the string frame metric then increases in the 
IR, As ^ (l/2)ln(2) at large z. As a result, at some 
position z* there is a minimum of Ag where the string 
equilibrates; see Fig. 2. This minimum is given by 

A;(z*) = 0, z*« 0.876 Ay^ (18) 

The simplest falling string example is a string extended 
in X direction and falling in z: X^ = {t,x,0,0, Z{t)). 
Plugging this into the action, one hnds the effective ac¬ 
tion of the problem for Z{t): 


dtdxe^^‘^^Wl - Z^ 


(19) 
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FIG. 2: The combination As(z) as a function of the holo¬ 
graphic coordinate « (solid) compared to its IR (large z) 
asymptotics (dashed). As(z) has a minimum corresponding 
to the equilibrium scale of the QCD string. 

which generates the equation of motion (EOM): 

d, . ( 20 ) 

Instead of solving this EOM, for one string one can use 
instead conservation of energy in our time-independent 
background. The Hamiltonian for this case is 

h = ( 21 ) 

\/l - 

By setting it equal to the energy E, one gets directly the 
first derivative of Z. 

As usual, motion in general occurs between two turning 
points in which the velocity vanishes Z = 0. When z = z*. 
and the energy corresponds to the minimum, the string 
simply “levitates” without motion at this point. The 
holographic image of this stationary string, calculated 
by standard rules, generates some stress tensor T^''(x) 
distribution, with Minkowski indices = 0..3, describ¬ 
ing a static QCD string. The integral / d^xT^°{x) per 
unit length in x is known as the QCD string tension 
Tg = (420MeE)^. Rather than predicting it, one can 
use it to fix the fundamental string tension Tf, [ 8 ]. 

Near the levitation point one can approximate the mo¬ 
tion by a harmonic oscillation with the frequency 

a;, = [2A"(z,)]VE (22) 

However, because the potential is quite asymmetric, the 
validity of such harmonic approximation is rather lim¬ 
ited. In general, oscillations toward large z - the IR ~ 
take longer time and reach higher values of z, as com¬ 
pared to motion in the UV direction. 

Let us now return to the discussion of the holographic 
image of the string on the boundary. Standard rules re¬ 
quire calculation of the bulk-to-boundary propagators for 
all bulk helds, which source the operators of the gauge 


theory. One of them - the metric 5 ^;^ - sources the stress 
tensor mentioned above. The characteristic scale 
of the string energy distribution in transverse plane is 
given by the mass of the corresponding mode, known 
as the tensor glueball M(2++) « 2 GeV. This (rather 
large) scale defines the scale of the QCD string width 
5y ~ 1/M(2++) « 0.1 fm. 

Elavorless scalars source the quark bilinear operator 
qq, related to modification of the quark condensate in the 
QCD vacuum by the string. The lowest scalar meson in 
QCD is much lighter than the tensor glueball M(0“*'''') = 
nia ~ 0.4 GeV, and thus it produce much wider image. 
This field will play prominent role in what follows. 

The main subject of the present paper is not the mo¬ 
tion of an individual string, but rather the interaction 
and collectivization of many. A convenient experimental 
setting is provided by pA collisions. In peripheral col¬ 
lisions, an incoming proton collides with one or a few 
nucleons, producing several strings. But in central colli¬ 
sions the number of interacting, or “wounded”, nucleons 
is rather large, ^ 20, producing at least Ng ~ 40 
strings. Transition from a dilute to dense regime can 
thus be experimentally studied as a function of the mul¬ 
tiplicity of the produced secondaries. 

Furthermore, in real heavy-ion collisions {pPb at LHC, 
dAu, He^Au at RHIC), the multiplicity of secondaries 
depends linearly on the rapidity, with the ratio between 
the positive and negative rapidity ends being of the or¬ 
der 2. The Lund string model explains this by the fact 
that only a fraction of the strings go from valence quarks 
to valence quarks in the fragmentation regions of light 
and heavy nuclei; all others end at the gluons located at 
intermediate rapidities. 

For geometric simplicity, however, we will ignore such 
complications, and consider ensemble of strings which are 
parallel and infinitely long, so that their motion can be 
reduced to 3-dimensional coordinate normal to the time 
and beam directions, referred to as the “spaghetti” state. 
The number of strings we use should be understood as 
a local number at mid-rapidity, which is where most of 
the collider detectors and measurements are done. Such 
approximations we have inherited from recent study [19], 
which first addressed the issue of string collectivization 
not in the context of AdS/QCD. 


III. THE SCALAR MESONS 

A. The Scalar Masses in the Model vs. 

Phenomenology 

Before we turn to our model calculations, let us for 
completeness present a brief summary of the relevant 
hadronic spectroscopy. The scalar flavor-singlet sector 
of hadronic spectroscopy we need to address is one of 
the most complicated ones; its detailed discussion can 
be found in publications of Particle Data Group and the 
vast literature on which it is based [24]. 
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The scalar meson channel includes one light state, 
known as the a meson, or /o(600) in the current PDG 
naming scheme. Multiple complex fits put its location 
around M + ir/2 « (400 - 600) + *(200 - 300) MeV, but 
with a large spread. It is followed by a group of relatively 
close states /o(1370),/o(1500), and /o(1710). In ge_neral, 
all of them are expected to be a mixture of uu + dd, ss, 
and the lowest glueball states. A number of arguments 
- narrow width, diffractive production, and absence in 
77 production channel - indicate that the middle state 
/o(1500) is mostly comprised of the glueball. Splitting of 
/o(1370) and /o(1710) is comparable to 2ms, and their 
decay channel suggests that the upper state is predomi¬ 
nantly a strange ss meson. 

The chiral symmetry breaking phenomenon is strongly 
related with the existence of the light scalar meson 
a — /o(600). The standard form of its interaction with 
quarks vacuum, aqq, shows that the vacuum expectation 
value (VEV) (a) is - up to a sign - nothing else but the 
“constituent quark mass.” Partial or full cancellation of 
this VEV implies local restoration of chiral symmetry. 

It is crucial to our project to include the back-reactions 
of the quarks on the gluonic observables, or use the 
Veneziano limit. The study of the spectrum of such 
models has been recently done in [14, 15, 21, 22]. In 
particular, we use the confined phase spectrum with zero 
quark mass and the holographic potential of class I, which 
was analyzed in [15]. To calculate the meson and glue¬ 
ball spectra in holography, we have to expand the action 
to quadratic order and solve the linear Sturm-Liouville 
problem for the fluctuation fields in the bulk space-time. 
The numerical integration of the excitation equations 
starts from the IR with normalizable boundary condi¬ 
tions, and then proceeds towards the UV boundary of 
space, where we impose normalizable boundary condi¬ 
tions. Mass values which cause “bad” (non-normalizable) 
solutions to vanish in the UV provide the spectrum of the 
model. The masses of the lowest five modes in the scalar 
channel can be read off of Fig. 3. The glueball spectrum 
has also been studied in the context of other holographic 
models, such as the Sakai-Sugimoto model, but in the 
limit of small number of flavors, [25]. 

As seen from this plot, and previously shown in Fig. 
7 (left) of [15], at physically relevant flavor parameter 
X = Nf/Nc « 1, the lowest non-flavored scalar is indeed 
significantly lighter than the others. Fig. 8 (left) of [15] 
puts its mass to about half of the p meson mass and in 
the phenomenologically expected ^ 400 MeV range. 

Next come a close pair of the second and third states, 
with mass ratios to the first one m^jmi « 2.6. Since 
in the calculation the strange quark is as light as u, d, 
there should not be a separate /o(1710) state, and this 
pair can be identified with a close pair /o(1370), /o(1500); 
at a; = Nf/Nc = 1 their splitting is also correct. Different 
x-dependence of the third state from others hints that it 
is indeed mostly a glueball, but this feature is not robust, 
as it depends on the details of the potential. 

The five lowest masses in units of the UV scale of the 


model are 


mi 

Auv 

m4 

Auv 


1.53, 

4.86, 


m 2 

ms 


3.54, 


ms 

Auv 


5.45. 


3.94, 

(23) 


Selecting the numerical value of Auv from the second 
and third state masses, which are narrow and therefore 
well mapped to phenomenology, we fix the absolute units 
in our model to be 


Ac/v = 387MeV, m,^ = 592MeV. (24) 

Then, of the masses of the higher states read 


m 2 = 1370 MeV, m 3 = 1525 MeV, 

m 4 = 1881 MeV, ms = 2019 MeV. (25) 


B. Mixing of Pure States 

In the language of QCD, there are two kind of unfla¬ 
vored hadrons: those made of glue - glueballs, and those 
made of quark-antiquark pairs - mesons (qq). In the 
string language these are, respectively, made from closed 
and open strings. In the AdS/QCD models those origi¬ 
nate from fluctuations of closed strings coming from the 
color branes or open strings from the flavor branes. 

In any of those pictures, states with identical quantum 
numbers can mix, and when the mixing is strong, the 
original designation of states loses its meaning and one 
can only follow evolution of those as a: = Nf/Nc changes. 
However, if the mixing remains relatively weak, such a 
distinction of states remains meaningful and can help one 
to understand what is happening, in particular in the 
small-magnitude of the string-string interactions that we 
focus on in this work. 

With such motivation, we have studied the mixing phe¬ 
nomenon in significant detail. The excitation equations 
were derived in [15], Eqs. (AlOO, AlOl), for the two 
gauge invariant combinations of scalars defined in (II A) . 
We use them in the following form, of two coupled equa¬ 
tions 


c" + HA) C + p(A) c' + z\A)^ + N, (A) (C - 0 = 0, 

(26) 

C” + HA) C + HA) H + HA) m^^ + N2{A) (^ — ^) = 0 , 

(27) 

where the coordinate used is the scale factor A instead of 
the usual AdS coordinate z and the derivatives are with 
respect to A. This choice of coordinates is that the nu¬ 
merical calculations close to the AdS boundary become 
substantially more accurate. The coefficients in equa¬ 
tions (26) and (27) are, expressed in the A coordinate: 
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z"(A) 


k{A) = ( k{z{A))z'{A) - ) , 


t{A)=t{z{A))z'{A)\ 



m 

Auv 


where the original coefficients (without the tildes) are 
given by the lengthy expressions in [15], Eqs. (A. 102 - 
A.107). We will not copy them here, but just comment 
that all parameters of those coefficients are fixed by the 
Lagrangian, and so, even with all the complexity of ex¬ 
pressions, everything is fully determined. Normalizable 
solutions, both in IR and UV, provide the scalar masses 
that we already discussed above. 

We defined a set of “zeroth order” states C,n '^,as 
the eigenfunctions of the Hamiltonian without mixing 


= 


0 \ 

0 H^J ■ 


(32) 


The eigenvalues are 
and 


Hc = - 


wc{A) f 

z'(A)2 




i{A) 


“unmixed masses squared” 


,(0)2 


(33) 

(34) 


The normalization weight for each is found by the stan¬ 
dard transformation to the Sturm-Liouville form of the 
uncoupled fluctuation equations, and all of those states 
are subsequently normalized to the unit norm. Hence, 
we have defined 


w^iA) = z'(A)2e/*(^) , w^{A) = . (35) 


In Fig. 3, we show the determinant of the UV 
boundary value of two linearly independent solutions of 
the scalar fluctuation equations as a function of mass 
squared. Zeros correspond to the normalizable solutions 
and denote the eigenvalues. The curve with closed circles 
corresponds to fully-coupled system, giving the “mixed 
mass squared,” while two other curves are for uncoupled 
equations. The lesson from this plot is that each mixed 
state is close in its mass to one of the unmixed states we 
use as a basis; this was our early indication that mix¬ 
ing effects are, in some sense, small. Looking at the Fig. 
3 more attentively, one finds an expected pattern of re¬ 
pulsive mass levels due to to mixing: close pairs of states 
move away from each other, the lowest state moves lower, 
etc. 


FIG. 3: The determinant of the UV boundary value of two 
linearly independent solutions of the scalar fluctuation equa¬ 
tions, versus the mass parameter in Auv units. The solid 
curve’s five zeros (red circles) indicate the normalizable solu¬ 
tions, and the corresponding masses are those of the lowest 
five mixed scalars. The two other curves (dashed and dot- 
dashed) correspond to unmixed equations as explained in the 
text. 




(Auv“’) 


(a) 



W((,zX„(zf 



(b) 


FIG. 4: The first three (a) meson and (b) glueball wave func¬ 
tions squared, multiplied by the appropriate weight. 
















To make it quantitative we proceed by defining the 
mixing part of the Hamiltonian 


V = 


_( 0 

0 


where 


^ wc(A) 

^ z'(Ay 

i{AY 


v, = - 


PiA)—-N,iA) 
q{A)^-N2{A) 


(36) 


(37) 


Using unmixed states as a basis, we define a vector 

(d°^) \ 

(d°^) / ’ 


vtl = 
^ nm 


where m,n = 1,2,3,... We then calculate matrix ele¬ 
ments Vnmn'm' = {'^nm\ V\^n'm’) of the mixing. Limit¬ 
ing ourselves to the subset of the first three glueball and 
meson states, we define the decoupled eigenenergies as 


nn' ^mm' 5 (38) 

where n,m = 1, 2,3 and 

and Enm is given by the uncoupled energies, correspond¬ 
ing to meson masses of (1.688, 3.709, 4.975) and glueball 
masses of (3.846, 5.370, 6.386) in units of A^y. The 
eigenfunctions of these first three meson and glueball 
states are shown in Fig. 4. 

The total Hamiltonian in the basis of takes then 
the (infinite) matrix form, of which the part related to 
first 3 states of each kind is given by the following 6x6 
matrix 


Hq + V = 

/ 17.6375 0. 0. -0.8702 0.2779 -0.1221 \ 

0. 34.3941 0. -0.4534 -1.4549 0.4314 

0. 0. 48.6373 0.2952 -0.1115 -2.4406 , . 

-9.2289 -2.5654 2.5144 3.3982 0. 0. ’ 

4.4342 -12.3857 -1.4334 0. 16.4069 0. 

\-0.5529 1.0396 -8.9644 0. 0. 29.5146/ 

in which the diagonal values are the unmixed mass 
squared. The eigenvalues of this matrix give the new 
(mixed) energy values, from lowest to highest, which are: 

Enew = (2.803,15.056,18.566,28.380,35.442,49.741), 

(40) 

which yields 

mnew = (1-533,3.553,3.945,4.878,5.451,6.458), (41) 


in units of Ajjy, where the first, second and fourth val¬ 
ues correspond to the three lowest meson states, and the 
third, fifth, and sixth correspond to the glueballs (similar 
to what is shown in Fig. 3). The coupled equations, on 
the other hand, gave the following masses for the first five 
mixed states: 1.533, 3.541, 3.943, 4.863, 5.447 in units of 
Ajjv 

The eigenvectors can of course also be readily found. 
In principle those can be used to predict specific decay 
modes of the physical states. We do not do this, but only 
draw attention to the lowest of these states, the mixed 
sigma. Its decomposition, for x = 1, over these first six 
unmixed states reads 

|cr) = .0587110) -b .0140 |20) - 0.0065 |30) -b 0.9932 |01) 

- 0.0075 |02) - 0.0015 |03) , (42) 

shown in Fig. 5 (a). One finds that the unmixed lowest 
meson state | 01 ) is very close to |cr), while the largest 
admixture to it is from the first glueball state |10). This 
statement however is not universal; while the second state 
vector is 

|cr') = 0.0888 110) - 0.0760 |20) - 0.0022 |30) - 0.0541101) 

- 0.9916 |02) -b 0.0074 |03) , (43) 

which is close to the second unmixed meson, only half of 
the third state is the unmixed first glueball 

|/o(1500)) = 0.5428 |10) -b 0.0599 |20) -b 0.0085 |30) 

- 0.339 |01) -b 0.7655 |02) -b 0.0287 |03) , 

(44) 

as observed in Fig. 5 (b) for x = 1. 

Note also that we have ignored many more states which 
can admix to these scalars, such as meson-glueball and 
two-meson states; accounting for them will reduce the 
“original” fraction of the glueball a bit more. 

The x-dependence of the mixing of the lowest mixed 
meson state and lowest mixed glueball state are shown 
in Figs. 5 (a) and (b), respectively. The plot depicts 
the squares of the mixing coefficients (found by taking 
the inner product of the mixed state with each unmixed 
eigenfunction and squaring). Particularly, Fig. 5 (b) 
explicitly shows that there is level-crossing at x larger 
than one, as was first noticed in [14]. 

The admixture of the glueball states to sigma is the key 
phenomenon we are interested in this paper. Indeed, the 
strings - as gluonic (dilatonic) objects - do not directly 
interact with the quark-related brane (the s scalar, in 
notations of [15]). One can express this fact with the 
string coupling vector [string) = ( 1 , 1 , 1 , 0 , 0 , 0 ) in the 
unmixed set of states. Its projection onto the physical 
sigma state is then 

(crlstring) « 0.066 (45) 

which provides the key small parameter of the problem. 
This was first extracted in [19] from lattice data of the 
two-string configuration. 
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FIG. 5: The square of the decomposition coefficients of the (a) 
lowest mixed meson state and (b) the lowest mixed glueball 
state. 


Sng = -Tf J dtdxe'^^‘’^^'>^x{z) ■ (48) 

Hence, at low velocities the NG Lagrangian, which 
sources the dilaton fluctuations is time independent. The 
source term in the C equation of motion, Eq. (26), is given 
by the y variation of the action, Eq. (48) 

-^et^o-^|^5(z-z')<5^(r), (49) 

where r = X 2 + X^. We now need to solve for the 
equations of motion sourced by the above delta-function¬ 
like perturbation. 


(Ho + V)aA,A') = -MA)fNGiA)SiA - H'), 

(50) 

where fNG{A) = ^ %\a)z^{A') ■ Primes on the functions 
denote the derivative with respect to H, and A' is the 
point where the source is located. This can be put in 
the matrix form described above, in which differential 
operators and delta functions are projected onto the set 
of unmixed states. We then have 



(nf+p^wc He \ 
He Hf+P^w^J 

^ ( -w^{A)fMG{A)5{A- A') 
\ 0 



(51) 


We expand the solutions and the source in terms of the 
C. Scalar Fields From a String Source ^( 0 ) eigenfunctions 


The Nambu-Goto action determines the string inter¬ 
action with the metric and the dilaton. The excitations 
are rewritten here 

dgMN = (JMN , (5$ = X. (46) 

The first-order correction of the Nambu-Goto action due 
to the above excitation fields is 


= E 


CnX 



fNG{A)6{A-A') = Y,di<\n. 


Sng = -Tf J 


Then, Eq. (51) becomes 



(G-^r^d^X^dpX^gMN) , 

(47) 

{niP ^ -b p^) -b cP J 

Hy’ 

1 

II 

J ^ 

■H-T 

(52) 

where we have defined Gap = g^.vdaX^^dpX''. The string 
is taken to be a X^ = {t,x,X 2 {t),X 3 {t),Z{t)). In the 

dO cPfmCn 

mP^+p^ 

JO (e^miCn) 

” mP^+p^’ 

(53) 


case of low string velocities, only the dilaton is sourced, 

and the first-order action is simply We define / = (C^l H^/4 ICn) leading to 
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In the non-relativistic limit, the equations of motion are 


c 


(C) 

t 


(OIVcIU (Cnlf^glCfc) \ 4^^ 

[rrSp'^ j + P'^ 


„(«) 


7(0 


Co' =- 




mP ^ 


(?)2 I 2 

ml 


(54) 


where dP = W|^(A')^„(A')/jvg(>1 '). The solutions then 
read 


Z{t) = -2dzAsiZ) - pzxiZ,X2,X3), 

X2{t) = --dx^x{Z, X 2 , X^), 

Mt) = -^dx,x{Z,X2,X^). (60) 


_ fNciA )W(;{A )Cri{A )Ck{A) 

' “ (C )2 n ^ 

mn’ +p2 

(Cnlnl 6 ) ( 6 inia) \ 

(mp^+p'^'^ (mp^+pp j 



(55) 


aAA') 


fNG{A')wc_{A')C,n{A')^l{A) (^^1 Vg ICn) 


00 2 


m 


(0 2 


(56) 


Yukawa potentials are given by the propagator, ^ 
exp(—m„r). Since the distances of interest are r ~ 1 fm, 
going from the first to the second scalar means an extra 
factor of ~ exp[—(m 2 — mi)r] ~ exp[—5] <C 1. Thus, one 
can safely ignore all scalar states except the first one. 

The solution in coordinate space reads 

az,x,z',x')=J2 J ^,e^^'^^-^'^'CiA{z),A{z'),p) . 

(57) 

Considering strings in the quasi-static limit, po = Oj ^-nd 
parallel in x^ direction with pi = 0 , one is left with 
only two transverse momenta P 2 ,P 3 over which to inte¬ 
grate. In the limit of p ^ mP, the first-order ^ field 
term gives a contact term in the x^x^ plane, so we keep 
only the second-order term to describe long-range inter¬ 
action. The potential is found to be a Yukawa potential 
in the x'^x^ plane of the form of the one used in [19] with 
additional dependence on the z-coordinate, which is gov¬ 
erned by the holographic dynamics. As shown above, the 
string couples to the y field. Selecting the gauge where 
•0 = 0, we have that x(z) = 


x(z,Y,r) = -(Ci|Yc| 6 ) ( 6 lnlCi)x 

fNG{z')wQ{z')l^l{z')C,l{z) A'{z) 


27r m\ 


(C)4 


$'(z) 


Ko{mP r). 
(58) 


We now consider string motion in the background po¬ 
tential and the x fluctuation. The string action is 

Sng = -Tf (^1 + ^x) ^1-Z2-A|-A|. 

(59) 


IV. CHIRAL SYMMETRY RESTORATION 


As was shown initially in [II] and further studied in 
[12, 13], chiral symmetry breaking is described by the 
dynamics of the bulk tachyon field, which is dual to 
the qq operator. The diverging IR asymptotics of the 
tachyon space-time lead to the breaking of chiral symme¬ 
try in the boundary field theory. The chiral condensate is 
read holographically by the UV asymptotics of the bulk 
tachyon field, Eq. (10). 

The Nambu-Goto string, in the gravity-dilaton- 
tachyon background, sources the scalar glueball fluctua¬ 
tion (x), which is also coupled to cr-meson (s), Eqs. (26, 
27). In the case of a large number of string sources, the 
excitations are strengthened and they can significantly 
change the background dilaton and tachyon fields. 

The fluctuation induced by one string is compared to 
the vacuum solution in Fig. 6(a), for a source string 
at the equilibrium point z*. As one can see, the z- 
dependence of both is similar, but the effect of one string 
is small compared to the background. 

Placing static strings around the equilibrium point, z*, 
we calculate the necessary density of strings that source 
an excitation field that cancels the vacuum solution close 
to the boundary. The density of strings is defined as 
the number of strings over the area that they occupy. 
We notice that since the string sources are static, the 
fluctuation field will also be time independent. The ratio 
of the two, for the critical string density ps = 10.85 fm“^, 
is shown in Fig. 6(b); one can see that it is close to unity 
for a large range of z. This density of coincidental strings 
thus cancels the vacuum solution, and consequently leads 
to restoration of the chiral symmetry, similar to other 
models [19, 26]; this restoration is unlike that of those 
studies, however, as they observe the modification of the 
chiral condensate within the flux tube. 

Since, in this study, we are concerned with the inter¬ 
action and dynamics of the strings, the strings will not 
always be at the equilibrium point of the potential. Due 
to the nature of the x fluctuation and its dependence 
on the source point z', the density of strings that are 
needed to source an excitation that cancels the vacuum 
field varies with location. This behavior is shown in Fig. 
7. One can see that if the strings are located further into 
the infrared (larger z), their individual contribution is 
larger, and thus, a lesser density of strings is necessary 
to restore chiral symmetry. 
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FIG. 6: (a) The (tachyon) excitation induced by a single 
string s is compared to the vacuum (background) field, (b) 
The ratio of the (tachyon) excitation produced by transverse 
density of strings pe = 10.85 fm“^ with the vacuum (back¬ 
ground) field. 


FIG. 8: The background potential, (a) without and with 
string-induced fluctuations, all placed at the minimum of the 
2 potential ( 2 *) with the denoted transverse density, and (b) 
induced by strings with density 11 fm“^, all placed at various 
points in the 2 coordinate (denoted 2 s). The r dependence of 
X is averaged out, leaving only the density dependence of the 
fluctuation. 


ps (fm 



0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 ' 

FIG. 7: The density of strings necessary to induce a (tachyon) 
excitation that cancels the vacuum (background) field as a 
function of placement in the 2 coordinate. 


The numbers given are based on a simple linear ex¬ 
trapolation from a single string. Obviously, when the ex¬ 
citation field that gets comparable to the vacuum field, 
the problem becomes nonlinear; the sources should be 


included in the vacuum equations of motion. To solve 
those is beyond the scope of the present work. 


V. COLLECTIVE EFFECT OF STRINGS IN 
THE BULK BACKGROUND 

In this section, we consider effects of strings, dis¬ 
tributed in the transverse 2-d plane with a certain den¬ 
sity, on the background potential. The string fluctuation 
deformation of the background due to different densities 
is shown in Fig. 8(a). We notice that the minimum of the 
potential no longer exists if a certain density of strings is 
present at z = 2 * (approximately 18 fm“^). 

Similarly, we also notice that, based on the form of 
x(z, z'), the further a string is placed from the holo¬ 
graphic boundary, the larger the induced fluctuation on 
the background. This behavior is shown in Fig. 8(b), 
in which the string fluctuation is sourced by strings with 
transverse density 11 fm“^. The strings are then placed 
at a certain point in the z coordinate, denoted Zg. We 
see that at Zg = 0.4 A^y, the fluctuation is almost negli- 
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gible - the background remains almost unchanged from 
the zero string case but if the strings are placed at 
Zg = 1.2 Ayy, the string induced string fluctuation is 
equivalent to a much larger density of strings placed at 
the original minimum of the potential; compare with Fig. 
8 (a). 

In the original background with no string sources, a 
string would be in a stable oscillating state when initially 
placed anywhere in the holographic direction. However, 
with the induced fluctuations of a number of strings, the 
string may be in, or move to, an unstable point of the po¬ 
tential and subsequently be pushed indefinitely far from 
the UV boundary (^ = 0). The configurations that ex¬ 
hibit this behavior are highly dependent on initial con¬ 
ditions: the number of strings, initial placement in the 
z-direction, and the transverse density. As can be seen in 
the equations of motion, the strings attract in the trans¬ 
verse plane, while they “fall” in the holographic coordi¬ 
nate. Therefore, when enough strings are close enough 
and induce a large enough fluctuation, on the order of 
magnitude of the original background potential, they col¬ 
lapse onto each other in the transverse plane and are 
repelled infinitely far from the UV boundary in the z- 
direction (i.e. “fall” into the IR). When this collapse 
occurs, the interaction of strings is no longer described 
by our model, as they collectivize into a black hole and 
back-reaction must be taken into account. 


VI. NUMERICAL STUDY OF MULTI-STRING 
MOTION 

To study the motion and interaction of a number of 
strings, we wrote a program in Python, utilizing its 
odeint differential equation module. The background 
functions, wavefunctions, and subsequent interaction 
functions were computed and analyzed in Mathematica 
and exported for use in the Python program. 

The initial conditions in the transverse (a:^, a;^) co¬ 
ordinates are the same as those given in [19]. In the z- 
direction, a simple Gaussian distribution was used with a 
standard deviation of 0.11 A[)^y. It is important to note 
that, in the z-direction, the interaction does not have 
to do with the actual distance between a fluctuation¬ 
sourcing string and any other string; the position of a 
sourcing string determines it’s effect on the background 
potential, which then changes the other strings’ behavior. 

Using these initial conditions, the number of 
“wounded” nucleons, iV^,, each of which creates 
two strings, was varied. The cases treated were 
= 2,3,4,5,6,7,8,9,10,11,13,15,18,20, and 25, 

each of which had the transverse distribu¬ 
tion described above and centered at z = 

0.44,0.55,0.66,0.77,0.88,0.99,1.1, and 1.2 A[)[^. Six¬ 
teen trials were run for each string number and initial z 
distribution. 

The timescale of interest is on the order of 1.5 fm/c. 
By the end of this time period, the strings break due to 



FIG. 9: Percentage of trials in which the strings collectivized 
before breaking (i.e. within 1.5 fm/c) as a function of number 
of strings and initial center of Gaussian distribution in z. 


Schwinger type qq production and subsequently disap¬ 
pear, as was mentioned above. 

The collectivization results are shown in Fig. 9. We 
see that the higher momentum states - those that are 
initially closer to the UV boundary - have a higher prob¬ 
ability to collectivize with a fewer number of strings in 
the time period of interest. Cases with A/j, = 20-25 col¬ 
lectivize regularly when initially close to the boundary 
but do not “fall” before breaking if they are near the 
minimum of the original potential. 

However, we also see that if they are far enough away 
from the boundary (around z = for 50 strings, 

for example), they once again begin to fall before break¬ 
ing due to their larger induced fluctuation on the back¬ 
ground potential (see Fig. 8(b)). 

An example of a configuration of strings that falls is 
shown in Fig. 10. This configuration of 30 strings cen¬ 
tered around 0.44 A/)^ collectivizes in approximately 1.4 
fm/c, which is when the strings all “fall” to the infrared 
(large z) region. We cannot see any behavior after 1.32 
fm/c in this particular trial because at least one of the 
strings falls beyond the IR cutoff of our numerics. The 
perturbation on the background caused by this string 
located at large z is large, since large z magnifies the ef¬ 
fect of each string on the background, as shown in Fig. 
8 (b). In this case, the background potential no longer 
has a minimum, so we see that as soon as one of the 
strings is far enough from the z = 0 boundary, the in¬ 
duced background fluctuation causes all other strings to 
subsequently follow the first to the IR. 


VII. SUMMARY 

The aim of this work was to study string dynamics in 
a holographic setting. We have chosen to use the V-QCD 






^3(fm) XjCfm) ^3(fm) 


13 



A'2(fm) 




t ★ 








3I-------13 

- 2-10 1 2 

^2 (fm) 


(a) 



(b) 



★ 

i 

★ 

★ 

% 

★ ★ 

*★ 

* ★ 
★ 

★ 


★ 






(c) 


FIG. 10: Snapshots in the transverse (left) and holographic (right) planes of a falling 30 string configuration initially centered 
at 2 = O.dlAyy at t = (a) 0 fm/c, (b) 1.0 fm/c, and (c) 1.32 fm/c. 
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holographic model, with background potentials of class I 
given in [15] . It is worth noting that we have not changed 
any parameters or potentials of the model in this work. 
In this sense, all we do is a straightforward calculation, 
in a framework fully defined by matching to other QCD 
features in previous works [13-15, 20-22], 

We studied the spectroscopy of the scalar hadrons by 
treating the full states as a mixture of uncoupled meson 
and glueball states, which are coupled perturbatively. We 
found that the perturbative spectra are in a particularly 
good agreement with the full numerical solution of the 
V-QCD model. We also found interesting level-crossing 
between glueballs and mesons. 

We then turned to QCD strings, the holograms of the 
fundamental strings in the bulk, levitating in the min¬ 
imum of the effective potential formed by gravitational 
and dilaton fields. The first qualitative result we obtain 
is that the QCD strings should have a “gluonic core,” 
associated with glueball mass, and a “sigma cloud,” as¬ 
sociated with the sigma meson mass. 

We then found the response of all fields due to the 
presence of a string by considering the coupling between 
the scalar glueballs and mesons to be perturbative. This 
approach is supported by the fact that the perturbative 
spectrum is very close to that of the fully coupled system. 
We found that a certain density of strings cancels the chi¬ 
ral condensate in the vacuum background, thus restoring 
the chiral symmetry. The magnitude of string-string in¬ 
teraction is related to the strength of the glueball-meson 


mixing. This is weak because, in the no-mixing approxi¬ 
mation, strings do not interact with quark-made mesons. 

We also found that corrections to the dilaton back¬ 
ground can modify it in such a way that the minimum of 
the effective string potential disappears. If so, there are 
no more QCD strings, since all fundamental strings are 
indefinitely falling into the IR direction. 

Finally, we considered motion of the multi-string sys¬ 
tem, initiated in a “spaghetti” configuration, in which 
they are all considered to be extended indefinitely in the 
“beam” direction, so that the motion is only happen¬ 
ing in 3 transverse coordinates: the plane transverse to 
the beam and the holographic coordinate z. Depending 
on the number of strings and their initial locations, one 
observes attraction between them, which can be strong 
enough to increase the local density by an order of mag¬ 
nitude or so, and therefore trigger the collective “falling” 
into the IR direction. This is physically interpreted as 
“string melting” into a QGP fireball, as is indeed ob¬ 
served in pA collisions of high enough multiplicity. 
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